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We find that the current radiative energy loss kernels obtained from the opacity expansion dra- 
matically violate the collinear approximation used in their derivation. By keeping only the lowest 
order in collinearity terms, models based on the opacity expansion have ^ 50% systematic uncer- 
tainty in the calculation of 7r° Raa in 0-5% most central RHIC collisions resulting in a systematic 
uncertainty of ~ 200% in the extracted medium density. Surprisingly, the inclusion of a thermal 
gluon mass on the order of the Debye screening scale affects Raa at only about the 5% level due to 
non-intuitive coherence effects. For some observables such as Raa, the effect of these uncertainties 
decreases with increasing jet energy; for others, such as the average number of radiated gluons, the 
effect is energy independent. We note that it is likely that the differences reported in the extracted 
values of medium parameters such as q by various jet energy loss models will fall within this collinear 
approximation systematic uncertainty; it is imperative for the quantitative extraction of medium 
parameters or the possible falsification of the hypothesis of weak coupling between the hard probes 
and soft modes of the quark gluon plasma medium that future radiative energy loss research push 
beyond the lowest order collinear approximation. 
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I. INTRODUCTION 

Jet quenching is a unique observable in ultra- 
relativist ic heavy ion collisions (URHIC) as high trans- 
verse momentum (high-p^) particles are the most con- 
trolled, calibrated, and direct probe of the fundamental 
soft degrees of freedom of the quark-gluon plasma (QGP) 
[l]. Rigorous falsification or confirmation of fundamen- 
tally different qualitative pictures of the basic physics of 
the QGP (e.g., whether it is strongly or weakly coupled 
0, its relevant degrees of freedom [MS], etc.) from com- 
paring theoretical predictions to high-pT data crucially 
requires a detailed understanding of both experimental 
and theoretical uncertainties. 

Once the qualitative picture is fixed, jet tomography 
[6l|7], the quantitative determination of bulk properties of 
the QGP through the study of the attenuation pattern of 
high momentum particles, becomes possible. Jet tomog- 
raphy requires both high precision experimental data and 
a precise theoretical understanding of partonic energy 
loss. Recent statistical analyses by PHENIX [8^ 9 as- 
suming infinite precision for theoretical calculations that 
assume a weakly-coupled QGP suggest that data are now 
certain enough to extract medium properties to within 
~ 20%. Given that no current perturbative energy loss 
model satisfactorily and simultaneously describes more 
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TABLE I: Comparison of the extracted dNg/dy for the work 
considered here compared to PHENIX 0-5% most central 7r° 
data. Note especially the nearly factor of 3 difference in the 
last three columns when using Rad Only (radiative energy 
loss only). Due to the use of the collinear approximation 
all three of these values should be considered equally valid 
determinations of the medium density. The assumed infinite 
precision for the elastic channel makes the convolved Rad+El 
extraction of dNg/dy suffer a smaller systematic uncertainty 
than from those studied in the radiative channel alone. The 
first two columns demonstrate the limited influence of the 
radiated gluon mass in measuring the medium density. See 
the text for more details. 



than one high-p^ observable (TOl HI] , such as the level of 
suppression and azimuthal anisotropy of the high-p^ glu- 
ons, light quarks, and heavy quarks as measured through 
the light meson p!2HT5] and non-photonic electron nu- 
clear modification factor [T6^-T9^, Raa{pt^4^)^ and the 
away-side suppression of jet-triggered hadrons, Iaa [20J, 
it is, perhaps, premature to claim that even a qualitative 
understanding of the jet quenching of the QGP exists. 
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Nevertheless, continuing with the assumption that the 
current perturbative quantum chromodynamics (pQCD) 
models accurately describe the physics, such quantitative 
extractions of medium properties can only be meaning- 
ful when the systematic uncertainties in the theoretical 
calculations and modeling are both known and included 
in the systematic uncertainties in the extracted quanti- 
ties. While there have been some previous qualitative 
estimates of the theoretical uncertainty stemming from 
the running of the strong coupling [21| i22j, the prob- 
ability leakage of the Poisson convolution [23l [24], and 
a phenomenological IR cutoff imposed to approximate 
the effects of a non-zero thermal gluon mass [25 , one 
of the main purposes of this paper is to begin the pro- 
cess of rigorously quantifying the combined experimental 
uncertainties and systematic theoretical uncertainties in 
extracted medium parameters. 

In particular, we investigate the theoretical uncertain- 
ties due to the collinear approximation and due to dif- 
ferent assumptions regarding the thermal mass of the 
radiated gluon in the opacity expansion approach [26- 
[30j for calculating the radiative energy loss of high-p^ 
partons in pQCD. We also asses how those uncertainties 
change when elastic energy loss is included [ST, 'SS^ . A 
summary of our results is shown in Table [l[ Most strik- 
ingly, we find that keeping only the lowest order term 
in collinearity results in ~ 200% systematic uncertainty 
on the value of the extracted gluon rapidity density for 
central Au-hAu collisions at top RHIC energies. While 
this paper does not quantify the uncertainties associated 
with the collinear approximation for other pQCD radia- 
tive energy loss models, [23, 26. 33, 34 , the uncertainties 
in those models are almost certainly similar to the results 
quoted here. In particular, the sizable discrepancies be- 
tween the extracted medium properties by different en- 
ergy loss groups [8 , 9 , 35 are probably within the current 
theoretical systematic uncertainty. 

There are, generally speaking, four main pQCD-based 
radiative energy loss formalisms applied to URHIC: opac- 
ity expansion (GLV) [26l Ell El EO] , multiple soft scat- 
tering (BDMPS-Z-ASW) [28",'36'-'45;, higher- twist (HT) 
j46ti49l. and thermal field theory (AMY) [50-53 . In this 
work we focus quantitatively on the opacity expansion 
approach; see Fig. [l] for a cartoon of the physics in- 
volved and visualizations of the key variables discussed in 
this text. Within the GLV formalism one derives an ex- 
pression for the single inclusive radiated gluon spectrum, 
dNg/dx, that is folded into a Poisson convolution [7 for 
the distribution of total radiated energy by the parent 
parton (see below). Like all radiative energy loss models, 
the formalism makes the eikonal approximation, namely 
that the parent parton has a sufficiently high energy that 
its path is approximately straight and that the interfer- 
ence from the away-side jet, of 0{1/E)^ may be safely 
ignored [54 . The formalism also neglects contributions 
from four-gluon vertices and assumes that the parent par- 
ton suffers independent, path-ordered collisions. These 
choices are justified in the opacity expansion approach by 
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FIG. 1: (Color online) Cartoon of the production, in- medium 
energy loss, and fragmentation processes that may occur per- 
turbatively for a high-pT parton produced in a heavy ion col- 
lision. The momentum labels are p for the outgoing parent 
parton, k for the medium-induced bremsstrahlung gluon, and 
q for the momentum transfer between an in-medium soft de- 
gree of freedom and the high-pT parent parton. Note the 
ordering of length scales displayed, <^ X <^ L. 



the assumption of a medium composed of Debye-screened 
scattering centers whose screening length <C A is 
much smaller than the mean free path of the parent par- 
ton [36 . Such an ordering of scales is consistent with 
thermal field theoretic estimates for /i and A [55 . 

Current evaluations of the diagrams resulting from the 
opacity expansion formalism drop a significant number of 
terms that are assumed to be small. These approxima- 
tions were made for analytic simplicity, but they do not 
seem to be inherently required. Specifically, GLV takes: 
(1) the collinear approximation, kr <C xE^ where k is the 
momentum of the radiated gluon and kr its component 
transverse to the motion of the parent parton; (2) the 
parent parton path length much longer than the gluon 
mean free path, L ^ A; and (3) the soft radiation limit, 
X <C 1, where x is the momentum fraction taken away by 
the radiated gluon (discussed further below). We note 
that all these assumptions are also made in the BDMPS 
and AMY formalisms^. HT makes a similar assumption 
to (2), does not make the soft gluon approximation (3), 
but does assume collinearity (1). After discussing (1) in 
greater detail we will briefly touch on (2) and (3) further 
below. 

This work focuses on the collinear approximation, 
kr <C xE, and, more generally, the effects of limiting 
the phase space into which gluon bremsstrahlung is al- 
lowed to radiate. The current derivations using the opac- 
ity expansion formalism yield a single inclusive gluon ra- 



^ In GLV and BDMPS (2) is used to neglect poles from propaga- 
tors multiplied by exp(— /xA^;) ~ exp(— /iA) <C 1, where Az is 
the distance between successive scattering centers; this approach 
is probably invalid for L < A ~ 1 fm. On the other hand AMY 
uses the central limit theorem in its Langevin approach and cor- 
rections are likely for L < 30A ~ 30 fm; this extra long path 
length is also required by the neglect of the interference between 
vacuum and in-medium induced radiation. 
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diation kernel, dNg/dxdkT, that knows nothing about 
the approximations used in its derivation; in calculat- 
ing dNg/dx the collinear approximation is enforced phe- 
nomenologically with a UV cutoff in the kr integration. 
We find that the dNg/dxdkr kernel maximally violates 
the assumption of collinear ity at small x. It is not surpris- 
ing, therefore, that dNg/dx^ and observables dependent 
on it such as Raa^ are highly sensitive to the 0{1) vari- 
ations of the cutoff one explores in estimating a system- 
atic theoretical uncertainty. Even worse, as we discuss 
further below, the use of only the lowest order collinear 
term means that collinearly equivalent definitions of x 
yield values of medium density extracted from data that 
differ by - 100%. 

Previous work [25 that investigated a phenomenolog- 
ical kT cutoff in the IR to approximate a thermal mass 
of the radiated gluon found a similarly strong sensitivity 
to the specifics of the cutoff. However, we find that with 
an explicit derivation of energy loss for non-zero gluon 
mass instead of a cutoff imposed by hand a posteriori^ 
surprising and non-trivial cancellations yield an energy 
loss that has little sensitivity to the exact value of rrig. 



II. OTHER SOURCES OF UNCERTAINTY 

Before returning to the primary focus of the 
paper — uncertainties introduced by the collinear 
approximation — it is worth enumerating here other po- 
tential sources of theoretical uncertainty that we do not 
attempt to quantify in this work but that, nonetheless, 
could be quite large. 

In this paper, we calculate radiative energy loss using 
the first order expression from the opacity expansion; i.e. 
the single inclusive gluon radiation spectrum is derived 
by scaling up diagrams with only a single in-medium scat- 
tering for the parent parton or its bremsstrahlung radi- 
ation, as is depicted in Fig. [l] by the average number of 
scatterings, L/X. The n*^ order in opacity explicitly eval- 
uates the interference terms, neglected in the lower order 
expressions, from diagrams with n in- medium scatter- 
ings. Within the set of approximations given above, and 
assuming the medium consists of Debye-screened static 
scattering centers, the GLV formalism yields a closed 
form expression for dNg/dx to all orders in opacity. Nu- 
merical study of dNg/dx [7l|27l|55] suggests that it does 
receive corrections from higher orders but that these are 
relatively small — around the ~ 30% level for relatively 
long paths of L = 5 fm and smaller for shorter paths — 
for RHIC- and LHC-like conditions. 

The Poisson convolution [71|36] of the dNg/dx kernel is 
an attempt to approximate the full probability distribu- 
tion of radiative energy loss, P(e), where the fraction of 
radiated energy is defined hy Ef = {l — e)Ei. The convo- 
lution assumes independent, incoherent emission of glu- 
ons; the effect of neglecting the interference between two 
or more emitted gluons is not currently known. There 
are also uncertainties associated with "probability leak- 



age" [71 [23] : in model implementations the parent parton 
energy is often not dynamically updated throughout the 
convolution leading to the possible violation of energy 
conservation (i.e. P(e) has support for e > 1). In this 
work, we assign the total weight from the convolution 
at e > 1 to a delta function centered at e = 1 (complete 
stopping), the so-called "non-re- weighted" approach [23 . 
Studies have shown that this procedure reproduces rea- 
sonably well distributions obtained from Poisson convolu- 
tion procedures that dynamically update the parent par- 
ton energy [56 . Moreover, for energy loss calculations 
based on the GLV formalism, the leakage tends to be 
small. When elastic energy loss is included this leakage 
is even smaller [57 . 

A realistic model of jet quenching in URHIC requires 
a weighted averaging over a large range of medium path 
lengths < L < 12 fm. QGP temperatures at RHIC en- 
ergies lead to a gluonic mean free path of A ~ 1— 2 fm [55] . 
Additionally, even in the most central collisions, a signif- 
icant portion of path lengths have L <2 fm, and the bias 
towards the surface from energy loss makes these shorter 
paths even more important [23] [32]. The assumption, 
L ^ A, is therefore violated for a large fraction of energy 
loss inducing processes. In the Raa calculations of this 
work, we simply apply the GLV formalism to all path 
lengths. There is an additional uncertainty associated 
with mapping the realistic medium density of heavy ion 
collisions — with its approximate Bjorken expansion and 
non-trivial, time-dependent transverse density profile — 
into the static, uniform "brick" problem in which the 
analytic formulae were derived. 

The small x approximation seems reasonable at RHIC 
energies where the single inclusive gluon radiation spec- 
trum dNg/dx peaks atx ~ jj./E < 0.05 <C 1. The small-x 
approximation should apply even better at LHC, where 
E will be larger by a factor ~ 10. However, the Poisson 
convolution widens the original dNg/dx spectrum — no 
matter how low the x value of the peak of the spectrum — 
thus distributing its weight out to larger values x. In this 
way the convolution introduces a sensitivity to the large 
X region of dNg/dx that is not well controlled due to the 
X <C 1 assumption. 

Probably the largest source of uncertainty, likely big- 
ger even than that due to the collinear approximation 
that we detail in this paper, will result from values cho- 
sen for and/or the running of as [2l][22j; since radiative 
energy loss varies as a^, it is highly sensitive to changes 
in the value of the coupling. One factor of as arises 
from the emission of the gluon bremsstrahlung while the 
other two result from interacting with the soft degrees of 
freedom in the medium; see Fig. [l] For these last two 
couplings it is not even clear what scale should be chosen 
for the running. Nevertheless, the momentum transfer q 
will always probe soft scales on the order of ja ~ 0.5 GeV 
> Aqcd ^ 0.2 GeV, and non-perturbative physics will 
be involved. The lack of a proof of factorization, whereby 
these non-perturbative contributions necessarily become 
small at high energies, is of course not particularly com- 
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forting. 



III. OPACITY EXPANSION ENERGY LOSS 
KERNEL AND x 

The first order in opacity result for tiie radiated gluon 
spectrum assuming massive quarks and gluons and using 
tiie Gyulassy-Wang model of static scattering centers [3] , 
to lowest order in x and kT/xE^ where = |k|, is given 
by [29. 



dNg _ Cna, L f d^q 2£k 
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Here, /3'^ = x^M'^ + m^, where M is the mass of the par- 
ent parton, and rrig is the mass of the radiated gluon; 
Cr = Cf for a parent quark, and Cr = Ca for a par- 
ent gluon. The corresponding result for massless par- 
ent quarks/gluons and/or massless radiated gluons is ob- 
tained by setting M and/or nig to zero in Eq. ([T]), respec- 
tively. After temporarily taking the qt = |q| medium ex- 
change momentum limit to infinity and making a change 
of variables (see [29 for more details) and additionally as- 
suming an exponentially decaying distribution, p{z)^ for 
the distance to the scattering center^, Eq. ^ becomes 



dNg SCRasjJ^L 
A 

krj-i ( krp 



dx 



7T X 
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dgrdkr 
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\3/2 ■ 



(2) 



We note that supposing the scattering center distribution 
to be uniform instead of exponentially decaying appears 
to make only a small difference in the gluon spectrum 

The integrand in Eq. (2| is both IR and UV safe. In 
principle, we could set the lower and upper qt and kr 
limits of integration to and oo, respectively. If the in- 
tegrand were exact, then it would have support only in 
the physical regions of the qr and kr integration space. 



In general it is necessary to specify the distribution in the dif- 
ference in distance between successive scattering centers. In the 
first order in opacity case, with only one scattering center, this 
difference in distance is between the scattering center and the 
production point. One can always take the production point to 
be at z = 0, and thus in this case the specified distribution is the 
distance to the first scattering center. 
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FIG. 2: Plot of dNg/dxdkr from Eq. ^ (solid curve) for 
a light quark with all masses set to 0, E = 10 GeV, L = 
5 fm, and representative values of /x ^ 0.46 GeV and A ~ 
1.25 fm for a medium density of dNg/dy = 1000 similar to 
RHIC conditions ^32j. Vertical lines depict the values of kr 
used as cutoffs to enforce coUinearity in Eq. Note that 
with X = 0.025 ^ fi/E^ dNg/dxdkr is large near kr ^ A;max, 
completely in contradiction with the collinear approximation. 
Dotted curve, from the unshifted integrand of Eq. ([T]), differs 
only slightly from the spectrum obtained from the shifted 
integrand of Eq. ^ (solid curve) . 



However, due to the small x and collinear approxima- 
tions, the integrand violates kinematic limits: one can 
clearly see that the integrand in Eq. Q has support over 
all qT and /ct (this is also true of the unshifted integrand, 
Eq. ([T])). We enforce physicality in the hope of better ap- 
proximating the exact result by restricting the qT and /ct 
integration region with cutoffs. A new result from this 
paper is a quantitative estimate of the systematic error 
that results from allowing 0{1) variations of these cutoffs 
on some observables and extracted medium parameters. 
Surprisingly, we find that dNg/dx is highly sensitive to 
variations in the /ct UV cutoff. 

For the qr integration we will make the usual choices 
<?min = and (7max = V^S/iE. The value for g^^ax is, in 
principle, the maximum momentum transfer allowed us- 
ing relativistic kinematics, though it has been pointed out 
that this is no longer correct with the change of variables 
q q+k used in the final steps leading to Eq. . Nev- 
ertheless, Eq. ([2| is rather insensitive to the qr cutoff [29] . 
and we show below explicit calculations of dNg/dxdkr 
using the shifted and unshifted integrands of Eqs. ^ 
and Q, respectively. These comparisons in Figs. [2] and 
|3] demonstrate that the shifted integrand with the un- 
shifted ^max well reproduces the results of the unshifted 
integrand, albeit with minor artifacts. 

Previous work by BDMS [25] found a strong sensitivity 



We thank Ulrich Heinz for finding this error in logic. 
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FIG. 3: Comparison of dNg/dxdkr spectra as a function 
of kr with x = 0.5 between GLV (all masses set to zero) 
and DGLV (non-zero radiated gluon mass) formulations us- 
ing both the shifted integrand of Eq. ^ (solid curves) and 
the unshifted integrand of Eq. ([T]) (dotted curves). The hori- 
zontal dashed line at is meant to guide the eye. The three 
vertical lines represent the choices of A;max described in the 
text. dNg/dxdkr better respects the collinear approximation 
for larger values of x as it has little weight near kr A^max- 
At small kr the radiated gluon mass, surprisingly, enhances 
radiation for DGLV due to coherence effects; at larger kr the 
mass has the expected effect of suppressing dNg/dxdkr- 



to a phenomenological IR cutoff in the kr integration of 
the BDMPS multiple soft scattering calculation, a resid- 
ual of the original vacuum radiation IR divergence. They 
also showed that the sensitivity of a quenching factor to 
variations in this cutoff decreased with increasing parent 
parton energy. The IR cutoff was imposed in order to 
approximate the influence of a non-zero thermal gluon 
mass; we find that including a gluon mass at the level 
of the gluon propagators [29 decreases the sensitivity of 
the energy loss and calculated Raa (see Table [l|) to the 
choice of gluon mass. This result is due to a non-trivial 
cancellation of effects that we will discuss further below. 

We do find, however, that the radiated gluon spectrum 
and resulting parton energy loss are extremely sensitive 
to the choice of the UV cutoff in the kr integration. Fur- 
thermore we find a sensitivity to the particular interpre- 
tation of X used in relating the components of k to the 
components of p. It turns out that the derivations of 
dNg/dx in the literature have used two different, albeit 
equal to lowest order in collinearity, definitions of x. We 
will use the light-cone normalization = p^ ±p^^ with 



mverse p 



,0,2 _ 



(p+ ±p~)/2, and denote p^ = E and 
= . With these conventions the two definitions of 
X are: (1) the fraction of plus momentum carried away by 
the gluon, x = x^ = k^/E^ [26l[27]; and (2) the fraction 
of energy carried away by the gluon, x = xe = k^/E |28j. 
In the usual notation where parentheses designate four- 
momenta, k = (/c^, /c^, k), square brackets light-cone mo- 
menta, k = k], and bold-faced variables are the 



transverse 2-vectors {kr = |k|), we find that the radiated 
massless gluon has on-shell momentum 



k = {xeE, y^{xEE)^-k^, k) = 



x^E+ 



,k]. 



(3) 

With these definitions of x+ and xe one may derive the 
exact relationships 



(4) 




Xe = x+ 



(5) 



where we have assumed that the initial parent parton 
momentum is P = {E^ E^ 0), and, hence, E~^ = 2E. 
From these formulas, it is easy to see that in the the 
collinear limit, x+ = xe^ so it is natural that the inte- 
grand in Eq. Q is the same when derived using these 
two different definitions of x: in one case terms of order 
kr/xEE are dropped whereas in the other terms of order 
kT/xj^E^ are dropped. Thus dNg/dx-^dkTdqT{x-^) = 
dNg/ dxEdkrdqT^XE) in the collinear limit. 

Two justifications for and derivations of a kT cutoff 
exist in the literature: (1) when interpreting x as xe^ 
^max = xeE keeps k always real (note that in the x = x+ 
representation k is real regardless of the value of kr) |43] : 
(2) when interpreting x as /c^ax = X-^E~^ enforces 
forward emission {k~^ > k~) [58]^. In this work, we will 
interpret /Cmax as enforcing consistency with the collinear 
approximation. It is, then, useful to determine that cut- 
off from a physical condition independent of the interpre- 
tation of x; we choose to set /Cmax by requiring the emit- 
ted radiation be within a cone of angle ^max centered on 
the direction of the parent parton. In this way we may 
control the collinearity of the radiation by varying the 
maximum angle of emission allowed. Simple trigonome- 
try relates 



^ ^ I tan(6>max/2), x 

y X£;£;sin(6>max), X 

In particular, when ^max = 7r/2 



Xe- 



Xj^E^ — 2x-^E, X — x+, 
XeE, x = xe, 



(6) 



(7) 



and we see that the two original justifications of the A^max 
cutoff are manifestations of the same physical condition. 



^ An even more restrictive cutoff for kr results if one also requires 
forward propagation of the parent parton [?]. Surprisingly this 
tighter cutoff, which forbids support for dN/dx for a; > 1 and 
therefore enforces energy conservation, leads to only a small 
change in dNg/dx |57| . and we will not use it here. 
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Name 


X 


Jacobian 


^max 




x+, rUg - /i/A/2 


x+ 


No 


7r/2 




x+, nig = 


x+ 


No 


7r/2 





X^{xe), ^max = 7t/2 


x+ 


Yes 


7r/2 





Xe, ^max — 7l/2 




No 


7r/2 





Xe, ^max = 7r/4 


Xe 


No 


7r/4 






radiating parton path length L and the coherence length, 



2xE 



'^coh 



TABLE II: Descriptions of the five main radiative energy loss 
implementations investigated in this work and their names 
when listed in diagrams. See the text for details. 



namely forward emission. It is worth emphasizing that 
with /Cmax given by Eq. ([7|), for any numerically equal 
value of x+ and x e the integrand in Eq. ([2| is integrated 
out twice as far in kr when interpreting x as x+ as op- 
posed to interpreting x as x^;. We note that a previous 
study [44 also varied the upper limit of the integra- 
tion. However that work investigated the radiation into 
jet cones of various sizes; here we are interested in quan- 
tifying the impact of regions of kr space over which there 
is currently little or no theoretical control. 

If the integrand of Eq. Q respected the assumptions 
that went into its derivation over the entire integration 
region, then it would have little weight for kT ~ xE. 
However, Fig.[2]shows that there are values of x for which 
dNg/dxdkr has large contributions from kr ~ xE. In 
fact, we should expect such a result: the Debye mass 
fi is the natural scale for /ct, and dNg/dxdkr will al- 
ways have a significant weight at /ct As a result, 
varying /cmax ~ xE will always lead to large changes in 
dNg/dx near x ^ /a/E. In the kinematic regime rele- 
vant for RHIC and LHC, numerical study suggests that 

(1) dNg/dx{x) reaches its maximum at x ~ jii/E and 

(2) dNg/dxdkrix ~ j^/E^ kr) reaches its maximum at 
kr ^ ^max; ioT X ^ fi/E the maximum value of dNg/dx 
rises quadratically with /cmax- These observations imply 
a dramatic sensitivity of dNg/dx — and any quantities de- 
rived from it — to the precise choice of /cmax- Increasing 
the radiating parton energy decreases the region of x over 
which dNg/dx is highly sensitive to the cutoff, x < ji/E. 
Naively, then, one might expect that increasing E would 
decrease the sensitivity of observables calculated from 
dNg/dx to variations of /Cmax- However, we will show 
below that even this expectation does not generally hold. 

Fig. [3] shows that, as expected, as x increases the 
assumption of collinearity, kr <C xE^ becomes a bet- 
ter approximation: dNg/dxdkT has little weight near 
kr ~ kmax^ and, therefore, for these larger values of x 
dNg/dx has less sensitivity to the exact choice of /Cmax- 
One may also observe in the figure the negative values of 
dNg/dxdkT for the massless case at small values of kT- 
These negative values are due to the destructive interfer- 
ence between the 0*^ order vacuum production radiation 
(the QCD analog of the usual beta-decay radiation spec- 
trum from QED) and the 1®* order medium-induced ra- 
diation. The interference is controlled by the ratio of the 



Qt 



(8) 



found by examining the cosine term in Eq. ([2| {qt en- 
ters Tcoh here due to the shift in integration variables 
Qt ^ Qt ^T in going from Eq. ([T]) to Eq. ([2|). For 
larger values of x the coherence length, Tcoh is longer, 
and the destructive interference more important; hence 
dNg/dxdkT is negative for small values of /ct for the 
massless results shown in Fig. [3] but not in Fig.|2j Includ- 
ing a non-zero mass for the radiated gluon reduces Tcoh 
and therefore the influence of the destructive interference. 
As a result, for small values of /ct the dNg/dxdkT spec- 
trum is actually enhanced for the massive as compared 
to massless case. At larger values of kT the non-zero 
thermal gluon mass exhibits the expected effect of sup- 
pressing bremsstrahlung radiation. Unfortunately, then, 
introducing a phenomenological /ct cutoff in the IR into 
a massless dNg/dx formula does not capture well the 
complicated dynamics resulting from allowing the radi- 
ated gluons to pick up a non-zero thermal mass from the 
medium. More readily visible in Fig. [3] than in Fig. [2] is 
the removal of the minor artifacts from the shift in inte- 
gration variable when evaluating the unshifted integrand 
of Eq. 0. 

The ambiguity in the interpretation of x has conse- 
quences beyond the differences in the corresponding lim- 
its of the kT integration. The gluon dNg/dx distribu- 
tions obtained from Eq. [2] using the two different x defi- 
nitions imply different gluon energy spectra and, thus, 
different amounts of jet quenching. These differences 
will clearly contribute to the systematic uncertainty in 
the interpretation of experimental data using the opacity 
expansion formalism and any other energy loss formal- 
ism that invokes the same or similar collinear approxi- 
mations. Furthermore, in order to better understand the 
consequences of the different physical pictures used in 
different formalisms — e.g. by comparing the gluon spec- 
tra produced by GLV (traditionally interpreting x as x+) 
and BDMPS-Z-ASW {x as xe) — one wants all results in 
terms of the same variables. As a first step we will ex- 
amine an "apples-to-apples" comparison between results 
derived from the two interpretations of x solely within 
the opacity expansion approach. 

There is an additional problem with energy loss calcu- 
lations that interpret x in Eq. (|2j) asx+: evaluation of the 
Poisson convolution P(e), where e is explicitly an energy 
fraction. Previous models that assumed the x+ interpre- 
tation in dNg/dx invoked the collinear approximation in 
order to directly use dNg/dx-^{x-^) unmodified as the in- 
put for the Poisson convolution. As we have shown this 
is a poor assumption. 

In order to compare the two interpretations of x 
within GLV and to quantify the effect of assuming 
dNg/dxj^ ^ dNg/dxE in finding P(e) we need to trans- 
form dNg/dx-^{x-^) to dNg/ dxEixE)^ which requires the 
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FIG. 4: Comparison of dNg/dxE{xE) for a 10 GeV quark 
using typical RHIC medium parameters and L = 5 fm (top), 
L = 1 fm (bottom) for the coUinearly equivalent expressions 
dNg/dx from Eq . (|2| with the x = xe interpretation and 
X-^{xe) from Eq. ([ofT^Also shown is dNg/dx from Eq. ([2| with 
X — xe and ^max = 7r/4 in order to explore the systematic 
uncertainty in the choice of /cmax- Inset shows results on a 
linear scale for clarity. 



FIG. 5: Comparison of dNg/dxE{xE) from GLV and DGLV 
(rUg = m/a/2) for a £^ = 10 GeV quark using typical RHIC 
medium parameters and L = 5 fm (top), L = 1 fm (bottom). 
Inset shows results on a linear scale for clarity. Note the 
enhancement of dNg/dx at L = 1 fm for the massive case. 



IV. QUANTITATIVE COMPARISONS 



use of a Jacobian: 



dNl 



-{xe) 



dxj^ 1 
dxE 2 



dxE dxj^dkT 



1+1 



f— )^ 

\xkEJ 



(9) 
(10) 



Note the change in the upper limit of integration in 
Eq. (|9| as dictated by the basic rules of calculus. The 
Jacobian, Eq. (10), is strictly greater than one and is 
singular as kr xeE; the competing effects between 
it and the numerically smaller /Cmax on dNg /dxE and 
energy loss will be discussed in further detail below. 



We now wish to quantify the effects of the different x 
interpretations, /Cmax values, and gluon mass treatments 
on dNg/dx and its derived quantities, such as P(e), the 
medium parameter dNg/dy extracted from Raa^ and the 
total radiated gluon multiplicity, (Ng). We will do so 
using five implementations of radiative energy loss; see 
Table im In the first three calculations we use the x = x+ 
interpretation and 6>max = 7r/2; in the first two we do not 
include the effect of the dx-^/dxE Jacobian. In the first 
calculation we assume radiated gluons acquire a thermal 
mass on the order of the Debye scale, rrig = ii/ \J2. The 
second calculation is the same as the first but with rrig 
0. The third is the same as the second but includes the 
Jacobian transformation from x+ to xe- In the last two 
calculations we adopt the x = xe interpretation from 
the start: in the fourth we take ^max = 7r/2 and the fifth 

6>max = 7r/4. 



8 



We show in Fig. [4]top (bottom) the dramatic diflFerence 
between the "equivalent" dNg/dxE spectra obtained us- 
ing the two interpretations of x and taking E = 1{) GeV, 
L = 5 fm (L = 1 fm), /i ^ 0.46 GeV and A ^ 1.26 fm 
as representative RHIC values. Also shown in the fig- 
ure is the large reduction in dNg/dx as 6>max is reduced 
from 7r/2 to 7r/4. That reduction in the spectrum demon- 
strates the importance of large angle gluon emission de- 
spite the assumption of collinearity. Fig. [s] (top) shows 
the dramatic reduction in dNg/dx with the introduction 
of a thermal gluon mass for the above parameter val- 
ues; the bottom panel in the same figure shows that for 
shorter path lengths, such as L = 1 fm, the inclusion of a 
thermal mass for the radiated gluon can actually enhance 
the emission of bremsstrahlung radiation. 

The only difference between the xe and Xj^{xe) 
curves, theoretically indistinguishable to lowest order in 



collinearity, in Fig. 
and higher. Fig. l6 



4 comes from terms of 0{{kT / xeEY 
quantifies the reduction of the effect 



of these terms as we decrease ^max from 90° to 30°. 
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FIG. 6: Ratios of gluon spectra obtained from Eq. 
dNg /dxE{xE), to that obtained from Eq. ^ with x = xe, 
dNg/dxE{xE), as a function of x^; for different choices of ^max- 
Smaller maximum opening angles produce a decreased sensi- 
tivity to the collinear approximation. 

It is also worthwhile to consider the "apples-to- 
oranges" comparison of the X-^{xe) curve in Fig. [4] and 
the rUg = curve in Fig. |5] (top). In particular, 
we observe that the former has a reduced peak height 
at X ~ fi/E but is larger for higher values of x. This 
generic behavior is due to a non-trivial interplay between 
the Jacobian, which always multiplies the dNg/dx-^dkr 
integrand by a number greater than 1, and the numeri- 
cally reduced /cmax = ^eE used in Eq. ([9|. Fig. [s] shows 
that the collinear approximation is better at larger values 
of X as demonstrated by the small values of dNg/dxdkr 
at kr ~ A^max compared to those at lower kr- The Ja- 
cobian, though, introduces an integrable divergence at 
kr = xeE^ so the X-^{xe) curve is systematically larger 
than the x+ curve at larger values of x. Since the total 



gluon yield, (Ng)^ must be unmodified by a change in in- 
tegration variables, the increased yield at larger x must 
come at the expense of gluon yield at smaller values of 
X] hence, the numerical reduction in k^ax wins in this 
small X region. This redistribution of probability means 
that (x) from X-^{xe) is always greater than {x) from x+. 
For any given set of medium parameters and path length, 
(Ng) must be the same for both the X-^{xe) and x+ im- 
plementations; since {x) is always larger for the former, 
using X-^{xe) will always produce a smaller Raa than 
For any given ^max, dNg/dxdkT with x = x+ is inte- 
grated out to a larger /cmax than for x = xe- Therefore, 
the x+ interpretation will always yield a smaller Raa 
than the xe interpretation. Using similar reasoning the 
Xe calculation with ^max = 7r/2 will always generate a 
smaller Raa than the xe with ^max = 7r/4 calculation. 
Based on these arguments we expect that the extracted 
medium parameter dNg/dy will be ordered from smallest 
to largest according to x^{xe), X-^, Xe with ^max — 

7r/2, 

Xe with ^max = 7r/4. The dNg/dy obtained via the sta- 
tistical analysis (described below) and listed in Table [l| 
demonstrate exactly this ordering. 

A crucial ingredient for any jet quenching calculation 
is P(e), the probability distribution for the fraction of en- 
ergy, e, radiated by a high-p^ parton. Described above, 
P(e) is obtained via a Poisson convolution of the sin- 
gle inclusive radiated gluon spectrum, dNg/dx^ assum- 
ing the independent emission of the multiple radiated 
gluons. P{e) can be decomposed into discrete and con- 
tinuous pieces. 



P{e) = P%{e) ^ P{e) ^ P^5{1 - e) 



(11) 



P^ is the probability of radiating no gluons (and hence 
no energy loss) and is given by 



eM-{Ng)). 



(12) 



P^ encodes the probability "leakage," the probability 
that the radiating parton loses a fraction of energy 
greater than unity, that results from the assumption of 
independent emissions used here. We show in Fig. [7| the 
P(e) generated from four of the dNg/dx curves investi- 
gated in this work (themselves displayed in Figs. [4] and 
|5|: x+ with rrig = 0, x+(x£;), xe with 6>max = 7r/2, 
and Xe with 6>max = 7r/4. The ordering of the aver- 
age fractional energy loss, (e), (represented by the ver- 
tical lines in Fig. [t]) is consistent with the qualitative 
arguments of the previous paragraph. As noted previ- 
ously {Ng) is identical for the x+ and X-^{xe) models; 
hence, P^ is also identical for the two. Previous evalua- 
tions of P(e) using the x+ interpretation approximated 
dNg/dxj^{xj^) ^ dNg/ dxE^XE)- We find a systematic, 
though small, difference between the (e) from this ap- 
proximation and that resulting from the correct transfor- 
mation of dNg/dx-^ to dNg /dxE] the modest variation 
in P(e) also produces only a ~ 10 — 20% change in the 
extracted dNg/dy. The continuous parts of P(e), P(e), 
have the expected ordering: from smallest to largest ac- 
cording iO Xe with ^max = 7r/4, Xe with ^max = 7r/2, 



This ordering follows simply from the integration of the 
dNg/dxdkT curves out to successively larger /cmax- The 
non-trivial ordering of P(e) for the x+ and X-^{xe) inter- 
pretations shows explicitly the redistribution of probabil- 
ity from smaller to larger e due to the combined effects 
of the Jacobian, Eq. (10), and the numerically reduced 
A^max of Eq. ([9|. 

Ultimately, we wish to quantify the effects on extracted 
medium parameters of: (1) the two collinearly equivalent 
interpretations of x, (2) varying /c^ax (via ^max), and (3) 
the thermal mass of the gluon. We will focus here on 
Raa{pt) measurements obtained from the 5% most 
central Au+Au colhsions at = 200 AGeV [8 . We 
follow WHDG [32 in the implementation of energy loss: 
leading order pQCD-based production spectra for glu- 
ons and light quarks, followed by in-medium energy loss, 
and, finally, KKP fragmentation into pions. Initial state 
effects such as Cronin enhancement and shadowing are 
neglected; as such the pr dependence of Raa here is 
stronger than in works that better describe the trend of 
the data [26 . 

So far in this paper we have only discussed medium- 
induced gluon radiation, and for one set of Raa calcu- 
lations we continue to use purely radiative energy loss. 
However, it is known that for pQCD in the kinematic 
regimes appropriate to both RHIC and LHC, elastic scat- 
tering leads to a significant contribution to the total in- 
medium energy loss Therefore, we will also 
calculate Raa using convolved inelastic and collisional 
loss. We use the Braaten-Thoma calculation [59l [60] 
of the mean elastic energy loss and use the fluctuation- 
dissipation theorem to estimate the width of the elastic 
energy loss distribution [32^ . We note that these assump- 
tions are a poor approximation for RHIC conditions be- 
cause the number of collisions is typically too small for 
the central limit theorem to apply. Hence, Gaussian dis- 
tributions do not do a good job of representing the ac- 
tual, highly skewed distributions involved; substantially 
improved results were derived in [56 . However, for the 
purposes of this paper we are not concerned with as- 
sessing the theoretical uncertainties resulting from the 
treatment of the elastic channel; this is a very interesting 
question in its own right. 

For the radiative energy loss calculation, we include the 
full multi-gluon fluctuations through the Poisson convo- 
lution. We also account for path-length fluctuations due 
to parent parton production points and trajectories in 
the medium for both inelastic and collisional loss with an 
approximate implementation of Bjorken expansion [32 . 
A standard Glauber modeling of the medium with a dif- 
fuse Woods-Saxon nuclear density function is used [32] . 
Specifically, hard production is assumed to scale with bi- 
nary collisions; the medium density is assumed to follow 
participant scaling. The strong coupling constant is held 
fixed at as = 0.3. We use thermal field theory to re- 
late /i = and A = l/pa, with a ex l/zi^ [32 . The 
only independent variable left is the single input param- 
eter dNg/dy^ the rapidity density of gluons in the pure 




0.001 



FIG. 7: Plots of P(e) for a 10 GeV quark in typical RHIC 
conditions with L = 5 fm (top) and L = 1 fm (bottom). P(e) 
obtained by convolving the single inclusive spectra dNg/dx 
shown in Figs. [4| and [sj Vertical lines represent the (e) for 
each P(e) distribution. Symbols represent the weight of the S 
functions at e = and e = 1. Inset shows results on a linear 
scale for clarity. 



glue QGP assumed here (increasing dNg/dy increases the 
medium density and decreases Raa)- The rigorous sta- 
tistical analysis of [8, 9^ was then used to determine the 
value and 1-a uncertainty of dNg/dy that yields a theo- 
retical curve that "best fits" the data given the experi- 
mental statistical and systematic errors. 

We show in Fig. [s] the PHENIX measurement of tt^ 
Raa{pt) for the 5% most central collisions and the best 
fit curves to the data for the five different implementa- 
tions of the radiative energy loss considered in this paper 
(see Table Il|) for both purely radiative energy loss (top) 
and convolved radiative and elastic loss (bottom). We 
note that, when compared at similar centrality bins, the 
STAR measurements of [6T| 7r++7r~ Raa systematically 
differ from the shown PHENIX tt^ Raa [8^ by as much as 
~ 50%. Such a discrepancy represents a potential source 
of additional systematic uncertainty in extracted medium 
properties that we do not attempt to evaluate here. 
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The results in Fig. |8] indicate that the extracted 
dNg/dy values for the massive {rrig = and mass- 

less {rrig = 0) emitted gluon cases are quite similar 
despite the large differences in dNg/dx for large path 
lengths, Fig. [s] (top). As we have shown this behav- 
ior is due to the interaction of the effects of a thermal 
gluon mass on the coherence length and the radiation 
kernel. As the path length decreases, the difference be- 
tween the massive and massless integrated dNg/dx de- 
creases until eventually the massive distribution exceeds 
the massless one. For central RHIC collisions the results 
shown in Fig. [s] (bottom) indicate this path length is 
~ 1 fm. As noted in [25] physical observables are quite 
sensitive to the specific numerical choice of the IR cut- 
off taken when used to approximate the inclusion of a 
thermal mass for the radiated gluon. The lack of sensi- 
tivity to the particular choice of gluon mass seen in this 
paper is a pleasant surprise. For the convolved energy 
loss, the extracted medium density is consistent for mas- 
sive and massless radiated gluons. For identical values of 
dNg/dy the massless case yields a smaller Raa] that the 
extracted medium parameter is consistent is likely due to 
the statistical analysis trading off goodness of fit in the 
normalization of Raa for slightly different pt dependen- 
cies. 

Addressing the primary goal of this paper, we ob- 
serve from the results in Fig. [S] that the extracted val- 
ues of dNg/dy for the radiative only energy loss models 
Xj^{xe)] xe, 6>max = 7r/2; and xe, 6>max = 7r/4 vary al- 
most by 200%. It cannot be overemphasized that the first 
two cases — for which dNg/dy varies by ~ 100% — are ex- 
actly equivalent under the collinear assumption used in 
the energy loss derivation. As discussed previously, for 
any particular ^max, Xj^{xe) produces the smallest Raa- 
Thus, Xj^{xe) with 6>max = 7r/2 represents a lower bound 
on Raa for fixed medium parameters. Unfortunately, 
it is not possible to rigorously define an upper bound 
to the theoretical uncertainty on Raa since decreasing 
^max makes Raa 1- We choose the values of 
observables calculated from the xe with ^max = 7r/4 im- 
plementation as a working definition of the upper bound. 
We consider this ^max as representing a reasonable 0(1) 
variation of the coefficient of the /ct cutoff that necessar- 
ily increases Raa for a given medium density; the sur- 
prise, of course, is just how sensitive Raa is to such a 
variation of the cutoff. It is worth noting that the larger 
dNg/dy values extracted using the radiative energy loss 
only models are difficult to reconcile with known RHIC 
dNch/dy multiplicities [62 . 

An important question that must be addressed is 
whether or not the sensitivity of extracted medium prop- 
erties to variations in the implementation of the collinear 
approximation decreases with increasing radiating parton 
energy. To address this question, we show in Fig.[9]a "toy 
model" quark partonic Raa{pt) calculation for L = 5 fm 
fixed path length and typical RHIC medium parameters 
used above. The sensitivity of the partonic Raa{pt) to 
different implementations of the collinear approximation 
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FIG. 8: The "best fit" curves of the five models discussed 
in the text to PHENIX measurements of the 5% most central 
Raa S '. top-including only radiative energy loss, bottom- 
including both radiative and collisional energy loss. Fits per- 
formed using the methods of [51 [9] . Uncorrelated errors are 
represented by bars whereas correlated errors are shown as 
grey boxes, and there is an overall scale uncertainty of ±12% 
not shown. 



decreases with increasing parton energy. For example, 
the :^ 100% difference at 5 GeV between the collinearly 
equivalent X-^{xe) and xe interpretations decreases to 
0^ 10% at 200 GeV. But, we observe that the number 
of emitted gluons, {Ng), remains sensitive to the uncer- 
tainties introduced by the collinear approximation; the 
~ factor of 2 difference in (Ng) shown in Fig. [9] is nearly 
energy independent. 

The convergence of the ratios of Raa{pt) with in- 
creasing energy can be understood at least partly as a 
consequence of the small x pileup in dNg/dx. For an en- 
ergy loss model with P(e) and a power law production 

spectrum dN/dpr cx p^*^^^^^ we have that Raa{pt) — 
/(I — e)^P{e;pT)de. Then for any two models, say 1 and 
2, 



^AA 



^((ei)(PT)-(e2)(PT)), (13) 



11 



and the ratio of modification factors automatically ap- 
proaches 1 as the mean energy loss decreases with in- 
creasing On the other hand, the constancy of {Ng) 
in pt is due to the surprising energy independence of the 
area under the peak in the dNg/dx distribution. As the 
energy increases, the distribution becomes more sharply 
peaked, and its normalization becomes more and more 
dominated by this energy-independent area. 
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FIG. 9: Ratios of a toy Raa model (top) and ratios of the av- 
erage number of emitted gluons, (A^), (bottom) as a function 
of parent quark energy for the five models discussed in the 
text. Collinear and gluon mass effects die out with increasing 
energy for the toy Raa but do not for {N). 



V. CONCLUSIONS 

In our introduction, we laid out the numerous approx- 
imations used in deriving the opacity expansion energy 
loss formulae. In one form or another, these assumptions 
are at the foundation of the four energy loss calculations 
that are frequently compared to RHIC data. In partic- 
ular, these calculations all rely on the lowest order term 
in a collinear expansion, where the small parameter is 
kT/xE. Unfortunately, we have found that the opacity 



expansion energy loss kernel dNg/dxdkT drastically vi- 
olates the collinear approximation for small values of x, 
the region of x most important for computing many ob- 
servables including the leading particle suppression and 
the average number of emitted bremsstrahlung gluons. 
While we did not explicitly check this violation for the 
other energy loss formalisms, it is highly likely they also 
will be similarly strongly affected by the large angle ra- 
diation that, by assumption of collinearity, is not under 
theoretical control. 

To leading order in collinearity, the two common defi- 
nitions of X used in energy loss calculations — x = x+, the 
fraction of light-cone plus momentum, and x = Xe^ the 
fraction of energy taken away by the radiated gluon — 
are equivalent. We found that for RHIC conditions, the 
"best fit" medium density dNg/dy extracted using these 
two different, but collinear ly equivalent, definitions of x 
varies by a factor of ~ 2. 

While one may use the X-^{xe) interpretation with 
^max = 7r/2 as a natural lower bound for the system- 
atic theoretical uncertainty of Raa, the upper bound 
is less obvious. We take the results from the xe and 
^max = 7r/4 as a working definition; in this case the ex- 
tracted medium density increases another ^ 50% over 
the Xe and ^max = 7r/2 model. We note that there is 
no sense of a "central value" or Gaussian distribution for 
this uncertainty band: the rigorous notion of a "best" in- 
terpretation of X or the "correct" UV cutoff for the lowest 
order collinear results does not exist prior to a calculation 
based on a more exact analytic derivation. 

The uncertainties we quote decrease significantly when 
elastic energy loss is included although the effects of un- 
certainty in the collisional channel were not considered 
here. While the leading particle Raa appears to suffer 
less systematic uncertainty at LHC energies, the uncer- 
tainty in other observables, such as the mean number of 
emitted gluons, are energy independent. Finally, the ef- 
fect of a thermal gluon mass on the extracted medium 
density is surprisingly small due to non-trivial coherence 
effects. This last result implies that the specifics of the 
short path length energy loss behavior are very impor- 
tant; future work should, therefore, go beyond the L ^ X 
approximation. 

One of the great debates over the past several years has 
been the so-called "discrepancy" of extracted medium 
parameters from the four energy loss models. An over- 
simplified description would be that the density found 
when comparing AMY, GLV, and HT to data is rea- 
sonably consistent while that found from BDMPS is a 
factor ~ 2 — 3 times larger [SI [9l |35] . It is natural and 
right to begin any calculation with assumptions about 
the relevant and irrelevant physics of the problem. How- 
ever, that the four models include and exclude vastly dif- 
ferent physics implies to us that any consistency found 
must be considered coincidental, even surprising. We see 
the path forward not in teasing out the origins of the 
extracted difference but in finding observables that can 
falsify the basic assumptions about the relevant physics 
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of the quark gluon plasma. Is the medium strongly or 
weakly coupled? What are its degrees of freedom? If en- 
ergy loss is perturbative, are parton interactions with the 
medium better approximated by many soft scatterings or 
a few hard ones? In order to answer these questions we 
must have not only well-controlled experiments but well- 
controlled theories. It is, thus, imperative that quanti- 
tative estimates be made of the systematic uncertainty 
introduced into theoretical results from the simplifying 
assumptions — distinct from the physics assumptions — 
made in the calculation. 

As discussed previously there is no notion of "one 
standard deviation" of theoretical uncertainty associated 
with the collinear approximation. Moreover, we have not 
quantified the consequences of the collinear approxima- 
tion for the other three energy loss models, although this 
is both an extremely interesting and clearly important 
problem. Nevertheless, we strongly suspect that should 
this uncertainty be quantified for the other models, then 
the different values of medium density so far extracted 
from data would be mutually consistent within the 
systematic theoretical uncertainties. As we have shown 
in this work, any hope for a quantitative extraction of 
medium density from high-p^ physics at RHIC using the 
GLV formalism requires a far more careful treatment 
of non-collinear radiation; this is almost certainly true 



for the AMY, BDMPS-Z-ASW, and HT formalisms, 
too. It is interesting that at LHC the leading particle 
suppression seems a rather collinearly safe observable 
whereas the average number of emitted gluons and the 
spectrum of soft gluons is not. This suggests that — at 
our current level of theoretical understanding — leading 
particle observables and jets measured using narrower 
cones may provide more sensitive tests of jet quenching 
at the LHC than, e.g., full jet reconstruction with large 
cones. 
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